Heisenberg frustrated magnets: a nonperturbative approach 
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Frustrated magnets are a notorious example where the usual perturbative methods are in conflict. 
Using a nonperturbative Wilson-like approach, we get a coherent picture of the physics of Heisenberg 
frustrated magnets everywhere between d — 2 and d — A. We recover all known perturbative results 
in a single framework and find the transition to be weakly first order in d = 3. We compute effective 
exponents in good agreement with numerical and experimental data. 
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Understanding the effect of competing interactions in 
three dimensional classical spin systems is one of the 
great challenges of condensed matter physics. However, 
after twenty five years of investigations, the nature of 
the universality class for the phase transition of the sim- 
plest frustrated model, the antiferromagnetic Heisenberg 
model on a triangular lattixie (AFHT model), is still a 
strongly debated question B. Due to frustration, the 
ground state of the AFHT model is given by a canted 
configuration ~ the famous 120nStructure - that implies 
a matrix-like order parameter o and thus, the possibil- 
ity of a new universality class. Experiments performed 
on materials supposed to belong to the AFHT universal- 
ity class display indeed exponents different from thq^e 
of the standard 0{N) universality class: for VCI2 
/3 = 0.20(2), 7 = 1.05(3), ZA= 0.62(5), for VBra 
a = 0.30(5), for CuFJID i: (3 = 0.22(2) and for 
Fe[S2CN(C2H5)2]2Cl ri: /3 = 0.24(1), 7 = 1.16(3). 
These results however call for several comments. First, 
the exponents violate the scaling relations, at least by 
two standard deviations. Second, they differ signifi- 
cantly from those obtained by Monte Carlo (MC) sim- 
ulations performed either directly on the AFHT model 
{f ~ 0.59(1), 7 ~ 1.17(2), /? ~ 0.29(1), a ~ 0.24(2)), and 
on models supposed to belong to the same universality 
class: AFHT with rigid constraints (i' = 0.504(10), 7 = 
1.074(29), /3 = 0.221(9), a = 0.488(30)), dihedral (i.e. 
1/3,2 Stiefel) models {1^ ~ 0.51(1)^7 ~ 1.13(2),/? ~ 
0.193(4), a ~ 0.47(3)). See Ref. § for a review, and 
references therein. Finally, the anomalous dimensions 
r] obtained by means of scaling relations is found to be 
negative in experiments as well as in MC simulations, a 
result forbidden by first principles for second order phase 
transitions O. All these results are hardly compatible 
with the assumption of universality. It has been pro- 
posed that the exponents are, in fact, effective exponents 
characterizing a very weakly first order transitktfv-fhe 
so-called "almost second order phase transition UllJ". 

From the theoretical point of view the situation is also 
very unsatisfactory since one does not have a coherent 
picture of the expected critical behaviour of the AFHT 
model between two and four dimensions. On the one 
hand, the weak coupling expansion performed on the 
suitable Landau-Ginzburg- Wilson (LGW) model in the 



vicinity of d = 4 leads to a first order niiase, transition 
due to the lack of a stable fixed point On the 

other hand, the low temperature expansion performed 
around two dimensions on the Non-Linear Sigma (NLct) 
model predicts a second order phase tiiansition of the 
standard 0(4)/0(3) universality class O. Since there 
is no indication that these perturbative results should 
fail in their respective domain of validity -i.e. for small 
e = A — d and small e = d — 2 - this situation raises two 
problems. First, and contrary to what happens in the 
non-frustrated case, one cannot safely predict the three 
dimensional behaviour from nai've extrapolations of the 
perturbative results. Although a direct compi|i4|ation in 



three dimensions, possible on the LGW model ^tJ'ta , can 
circumvent this difhculty, such an approach misses a sec- 
ond fondamental problem: the incompatibility between 
the symmetries of the NLcr and LGW models. Indeed, 
the renormalization group flow drives the NLcr model ac- 
tion towards an 0(4) symmetric regime, more symmet- 
ric than the microscopical system, a phenomenon that 
cannot occur within all previous treatments of the LGW 
model (see ref. [|l^ and below) . The LGW model is there- 
fore unable to find the 0(4) behaviour which Jias been 
nevertheless observed numerically in d = 2 c3. This 
raises serious doubts on the perturbative analysis of the 
LGW model away from d = 4. Reciprocally, the pertur- 
bative analysis of the NLcr model, based on a Goldstone 
mode expansion, predicts an 0(4)/0(3) fixed point ev- 
erywhere between d — 2 and d = 4, as for the = 4 
ferromagnetic model, in contradiction with the pertur- 
bative LGW results and the experimental and numerical 
situation in d = 3. All this suggests that non perturba- 
tive features could play a major role and thus imposes to 
go beyond the standard perturbative approaches. 

In this letter we realize this program hv using the Wil- 
son renormalization group framework Ell. We obtain a 
coherent picture of the physics of the AFHT model ev- 
erywhere between d = 2 and d = 4. We find that the 
fixed point expected from the NLcr model approach ex- 
ists indeed in the vicinity of d = 2 but disappears below 
- and close to - three dimensions. The transition for 
AFHT in d = 3 is thus weakly first order contrary to t 
different predictions of both a new universcdijy class 
and an 0(4)/0(3) second order behaviour O. We get 
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effective exponents compatible with the numerical and 
experimental data quoted above. For generalization to 
N > 4-component spins, we find the transition in d = 3 
to be second order with exponents in good agreement 
with recent extensive MC simulations - contrary to those 
found from three loop Padc-Borcl rcsummed series O . 

Our afipimach relics on the concept of effective average 
action , Ffc [(j)] , which is a coarse grained free energy 
where only fluctuations with momenta q > k have been 
integrated out. The field (j) corresponds to an average or- 
der parameter at scale fc, the analog of a magnetization 
at this scale. At the scale of the inverse lattice spacing 
A, Ffe^A is the continuum limit of the lattice hamilto- 
nian obtained, for example, by means of an Hubbard- 
Stratonovich transformation. On the other hand, the 
usual free energy F, generating one particle-irreducible 
correlation functions, is recovered in the limit k ^ 0. 
The fc-dependenrjp .of F^ is controlled by an exact evolu- 



tion equation 



dTk 
dt 



dt 



(1) 



where t = In/c/A. The trace has to be understood as 
a momenta integral as well as a summation over inter- 
nal indices. In Eq.(^, Rk is the effective infrared cut- 
off which suppresses the propagation of modes with-jno- 
menta q < k. A convenient cut-off is provided by EJ^E£l: 
Rk{q) — Zq'^ /{exp{q'^ /k'^) — l), where Z is the field renor- 
malization. In Eq.(^, F^^^ is the exact field- dependent 
inverse propagator - i.e. the second derivative of F^,. 

The effective average action Fj. is a functional invari- 
ant under the symmetry group of the system and thus 
depends on all the invariants built from the average or- 
der parameter. In our case, it is well known that the or- 
der parameter is a set of two vectors (pi and 02 that can 
be gathered in a real N x2 matrix (pab for A^-component 
spins o . The symmetry of the system is the usual spatial 
rotation group 0{N) times a 0(2) corresponding to the 
symmetry of the underlying triangular lattice EJ . This 
0(2) is realized on (pab as a right 0(2) "rotation" that 
turns the (pi into each other. There are two independent 
0{N) (g) 0(2) invariants built out of (pab- /? = Tr *(f>p and 



iTr(* 



The exact effective average action involves all the pow- 
ers of yO, r and of derivative terms, and so Eq.(|l|) is a 
nonlinear functional equation, too difficult to be solved 
exactly in general. We therefore need to truncate it. 
One possibility is to keep in F^ only the momentum 
(i.e. derivative)-independent part, an approximation 
called the Local Potential Approximation (LPA). In the 
case of frustrated magnets, this has been considered by 
Zumbach llJTij . This approximation however misses the 
field-renormalization and worse, as described below, the 
phenomenon of enlarged symmetjra around d = 2 found 
perturbatively in the NLcr model til . This does not mean 
that this approximation is not useful: it is simply, in 
essence, unable to answer the question of the matching of 



the different perturbative approaches. Another trunca- 
tion is however possible which preserves this possibility: 
it consists in an expansion of F^ around its minimum in 
order to keep a finite number of monomials in the invari- 
ants p and T while including the derivative terms allowing 
to recover the different perturbative results. We choose 
the simplest such truncation: 
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where {w, A, k, ^, Z} are the coupling constants which 
parametrize the model. All terms but one - the "current 
term" {eab<pcaS^ <pcb)'^ - are very natural and correspond 
to those appearing in the usual LGW action that real- 
izes the symmetry breaking scheme of frustrated mag- 
nets. Indeed for A and /i > 0, the minimum of the action 
is realized by a configuration of the form 0^*" = \/K.Rab, 
where Rab is a matrix built with two orthonormal A^- 
component vectors. The symmetry of this minimum 
is a product of a diagonal 0(2) group and a residual 
0{N — 2) group. The symmetry breakisg scheme is thus 
0{N) (g) 0(2) ^ 0{N - 2) O 0{2)diag Note that for 
0afc — (t>'ab^^ one has: p = 2k and t = so that Eq.(||) 
corresponds indeed to a quartic expansion around the 
minimum. The spectrum in the low temperature phase 
consists in 2A^ — 3 Goldstone modes and three massive 
modes: one singlet of mass mi = kA and one doublet of 
mass 7712 = which correspond to fluctuations of the 
relative angle and of the norms of the two vectors <pi and 

4>2- 

Without the current term, the truncation Eq.(|^) is 
however not sufficient in our case. This term plays a cru- 
cial role since, for A^ = 3, it allows the model to enlarge 
its symmetry from 0(3) ® 0(2) to 0(3) ® 0(3) - 0(4) at 
the fixed point around|-d.= 2, leading to the well known 
0(4)/0(3) behaviour O. The current term is system- 
atically discarded in the perturbative treatment of the 
LGW model around four dimensions, for the - correct - 
reason that it is power-counting irrelevant. Here we can 
include it in our ansatz since it is anyway present in the 
full effective action F^ and, in fact, must include it since 
it becomes relevant somewhere between two and four di- 
mensions. The formalism we use is in charge to decide 
where it is important. 

Let us emphasize that the effective average ac- 
tion method leads to non trivial and/or new results 
even within a quartic truncation of F^. One cap 
mention the Kosteditz-Thouless phase transition fEII, 
low energy QCD £i-, the abelian Higgs model and 
superconductivity Ol3', etc. The accuracy of the re- 
sults thus obtained depends on two main features: i) the 
smallness of the anomalous dimension 77 and ii) the fact 
that the thermodynamics of the system is controlled by a 
unique minimum of Ffc. Note finally that this technique 
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has been successfully employed in the case of the princi- 
pal chiral model to solve a conflict between perturbative 
approaches B , similar to what is studied here. However 
we stress that in the principal chiral case, there was no 
conflict between the symmetries of the LGW and NLcr 
models. 

The flow equations for the different coupling constants 
K, A, fi, u! and Z are derived by using Eq.(|l|) and Eq.(||) 
along the same lines as in [ p8[ . The explicit recursion 
equations are too long to display and not particularly il- 
luminating (see [^). Moreover, they require a numerical 
analysis, apart in d = 2 + e and in d = 4 — e where, as we 
now see, they get analytically tractable. 

The physics around two dimensions. Around two di- 
mensions, one expects that the perturbative "Goldstone 
mode" expansion of the NLcr model works well. In the 
Goldstone regime, the fluctuations of the modulus of 
and 02 and of their relative angle are frozen. This cor- 
responds to the large mass limit mi^, m2r — > oo. In this 
limit, our equations greatly simplify since the coupling 
constants divide in two sets {k, w, Z} and {A, /i} that do 
not mix. We only quote here the flow equations for the 
renormalized coupling constants of the flrst set: 



dKr . , „ ^ N -2 1 

= -(d _ 2 + rj)nr + + -7— r 

at ZTT 47r(l -f- KrLOr) 

— = (-2 + d + 2r])ujr+ 

1 + Kr-LUr + (N - 1)k'^UJ^ + (N - 2)k3cj3 
2K^7r(l + KrUJr) 

d\nZ 3 + iKrUJr + 2kIujI 



dt 



4Kj.7r(l -|- KrOJr) 



(3) 



These equations admit a fixed point for any > 2 of 
coordinates Kr — uJr — while Xr, ^J■r — est. The 
masses m*^, 7712^ are thus very large, proving the consis- 
tency of the limit. In fact, modulo the change of vari- 
ables: 771 = Kr and 'i]2 — 2Kr(l + KrUJr) the cquatious for 
Kr and uJr E^re exactly those obtained at pne-loop in the 
perturbative analysis of the NLcr model O. For N ~ 
they admit a fixed point for which the model is 0(4)- 
symmetric. 

Let us now recall how this phenomenon of enlarged 
symmetry for = 3 can be understood directly on the 
partition function. At the fixed point, the potential gets 
infinitely deep so that one recovers the hard constraints 
of the NLcr model: 4>i _L 4>2^ and 0f = 02 = ^or 
N = this allows us to rewrite the current term as the 
kinetic term of a third vector, the cross product of the 

two others: {eah(l>ca^ 4>cbY « (V^s)^ with ^3 = 01 A 02- 
The order parameter of the system is then a trihedral of 
orthogonal vectors (0i,02,03)- Thus contrary to what 
could be expected from a nai've expansion in powers of the 
fields, the current term plays a role as important as the 
usual kinetic terms. At the fixed point, ujr takes a value 



such that the three vectors play a symmetric role and 
the symmetry breaking scheme is 0(3) ® 0(3)/0(3) ~ 
0(4)/0(3) instead of 0(3) ® 0(2)/OC2^uch a result 
is of course missed within the LPA Therefore, 
the presence of the current term docs not only improve 
the accuracy of the calculation, it is necessary for its 
consistency. 

The physics around four dimensions. Around four di- 
mensions, we have expanded our equations at leading 
order in the coupling constants Xr and fir- At this or- 
der the current term decouples and we are left with the 
following equations for the quartic coupling constants: 



dXr 

~df^' 



4 + d)Xr + -^{AXrfir + 4/^^ + A^ (4 + A^)) 



^ = (-4 + 



167r 

+ ^{&Xrflr + Nfll). 



(4) 



They are i^^ose obtained at one loop from the LGW 
approach EJ . These flow equations admit a stable flxed 
point for A^ > A^^, ~ 21.8, attesting that the phase tran- 
sition is second order. For N < Nc the transition is first 
order since no fixed point is found. 

To higher orders, Nc depends on the dimension. In 
d = 3, three loop calculations ptpsummed d la Pade- 
Borel predict N^d = 3) = 3.91 H . Note however that 
this calculation exhibits unusual behaviours compared to 
the 0{N) case: the coefficients of the series do not de- 
crease mOjtiotonically and the series themselves are not 
alternate 113 . These features reveal the poor summabil- 
ity of the series. Finally, in the A^ = 6 case, for which 
the transition is second order, the predictions based on 
a Padc-Borcl rcsuxnmation, which provides = 0.575 
and 7 = 1.121 O, are in clear disagreement with re- 
cent numericalrsimulations, for which v = 0.700(11) and 
7= 1.383(36) y. 

From this point of view our approach has several avan- 
tages: first, since it matches with the one loop perturba- 
tive results in d = 2 and d = 4 it is likely that the error 
does not vary much with the dimension - a fact that 
has been confirmed in the 0{N) case for which the pre- 
cision for a given truncation is almost uniform with d. 
Second, it does not rely on a Pade-Borel resummation 
and therefore is free of the above mentionned problems 
of convergence. Of course, our results will change while 
improving the ansatz Eq.(^) by incorporating terms of 
higher order in fields and derivatives. However, all cases 
already treated within the average action method suggest 
that the lowest order approximation gives fairly good re- 
sults, even with this crude approximation. For example, 
in the ferromagnetic 0(3) model, one finds v = 0.703 c3 
which has to be compared to the six loop resummed 
perturbatioa, series in three dimensions which provide 
ly = 0.705 H. 

The physics between two and four dimensions. Let us 
first study the fate of the fixed point found analytically in 
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d = 2 + e for A'' = 3. By numerically integrating the flow 
equations, we find that this stable 0(4)/0(3) fixed point 
describes a smooth trajectory in the coupling constant 
space while d is increased. Our fiow equations actually 
admit another - but unstable - fixed point, which moves 
toward the stable fixed point as the dimension is in- 
creased. At a critical dimension dc ~ 2.87, the two fixed 
points collapse and disappear. Above dc, no other stable 
fixed point is found and we conclude that the transition 
is first order in d = 3. We thus show that the 0(4)/0(3) 
fixed point obtained from the NLtr model plays no role in 
the three dimensional physics of frustrated magnets, as 
conjectured for example by Jolicoeur and David M and 
Dobry and Diep & . We also discard the possibility of a 
new universality class conjectured on the basis of a naive 
extrapolation of the e = 4 — d calculation Gl'Bl . The prox- 
imity of dc with d = 3 however let open the possibility 
of a very weakly first order phase transition with effec- 
tive critical exponents. This behaviour manifests itself 
in our equations by the existence of a minimum around 
which the RG fiow slows down. This characterizes a very 
large, although finite correlation length ^. A rough es- 
timate of this correlation length - a few hundred lattice 
spacings - indicates that a pseudo-scaling behaviour can 
be observed although ^ is not large enough to ensure a 
true universality. This could explain the broad spectrum 
of effective critical exponents found in experiments and 
numerical simulations. Although the fiow equations do 
not have a fixed point, we are able to compute effec- 
tive exponents by linearizing the flow equations around 
the minimum. We recover here the phenomenon of "al- 
most secoBd pjjder phase transition" first introduced by 
Zumbach tJliJ within the LPA. To get accurate results 
we have to take into account the 0^-like terms in our 
ansatz. We find: ly = 0.53, 7 = 1.03 and /3 = 0.28, 
which lie in between the various sets of exponents found 
experimentally and numerically (see above). jJ'tp, com- 
parison Zumbach found j/ ~ 0.63 in the LPA Ollj, the 
difference being mainly due to the anomalous dimension. 

Finally, we find a true fixed point in d = 3 for N larger 
than a critical value Nc{d = 3) ~ 4. For iV = 6, we 
get V — 0.74 and 7 = 1.45 which compare well with the 
Monte Carlo data = 0.700(11) and 7 = 1.383(36) '1^ . 
They are close to the LPA results, where v = 0.76 O, 
and much better tkan those obtained by a three-loop cal- 
culation in d = 3 IIj (see above) . We have checked that 
our exponents do not vary significantly when monomials 
of order six in the fields are included in the ansate Eq.(^). 

To conclude, using a non perturbativc method, we 
have reached a global understanding of frustrated Heisen- 
berg magnets including a matching between previous per- 
turbativc predictions and a good agreement with ex- 
perimental and numerical data. It remains to under- 
stand the very origin of the disappearance of the NLtr 
model fixed point. The role of non trivial topological 
configurations can be invoked. One can hope a com- 
plete understanding of this point through the average ac- 
tion method which successfully describes the Kosterlitz- 



Thouless phase transition EJ . 
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